Diagnostics of plasma photoemission at strong coupling 
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We compute the spectrum of photons emitted by the finite-temperature large- A'^ SU{N) J\f — 
4 supersymmetric Yang-Mills plasma coupled to electromagnetism, at strong yet finite 't Hooft 
coupling. We work in the holographic dual description, extended by the inclusion of the full set 
of 0{a^) type IIB string theory operators that correct the minimal supergravity action. We find 
that, as the t' Hooft coupling decreases, the peak of the spectrum increases, and the momentum of 
maximal emission shifts towards the infra-red, as expected from weak-coupling computations. The 
total number of emitted photons also increases as the 't Hooft coupling weakens. 



Introduction. The analysis of data from heavy ion 
colhsion experiments at RHIC and LHC indicates that 
the quark-gluon plasmas (QGPs) produced are in the 
strongly-coupled regime [Ij , where the 't Hooft coupling 
governing the interactions of the microscopic constituents 
of the plasma is larger than one. Being electrically- 
charged, these microscopic constituents will emit pho- 
tons. The number of photons emitted with a given mo- 
mentum, i.e. the photoemission spectrum, yields valu- 
able information about the structure of the plasma. A 
theoretical study of this spectrum at strong coupling is 
therefore an essential step for investigating the QGP. 
Gauge/string duality lends itself perfectly to such a com- 
putation, because it allows the investigation of strongly- 
coupled gauge theories in terms of their weakly-coupled 
supergravity dual description Although there is no 
complete string theory dual model which accounts for all 
the relevant properties of QCD, the microscopic theory 
governing the behavior of plasma produced at RHIC and 
LHC, one can attempt to approach the real world using 
the holographic dual of the large— SU {N) A/" = 4 super 
Yang-Mills (SYM) plasma. Moreover, holography stipu- 
lates that A in the gauge theory maps to a'~^ in the string 
dual, so that the minimal supergravity {i.e. the zeroth- 
order string theory description in a small-curvature ex- 
pansion), obtained for a' — >■ 0, corresponds to the gauge 
theory at A — >■ oo. With these two caveats in mind, 
the work of Q considered two limits: photoemission 
from infinitely strongly-coupled SYM plasma, tackled via 
minimal type HB supergravity on the AdS-Schwarzschild 
black hole AcISbh x , and photoemission from weakly- 
coupled SYM plasma, computed using the traditional 
tools of perturbative quantum field theory. The real- 
world QGP lies somewhere in between these two illumi- 
nating yet unrealistic regimes. Our aim in this letter is to 
compute the photoemission rate of A/" = 4 SYM plasma 
at large finite A. We work in the holographic dual ex- 



tended by the inclusion of the full 0{a'^) type HB string 
theory corrections to the supergravity action. The rela- 
tion A ^ a'^"^ immediately dictates that the corrections 
to the X ^ oo result start at 0{\~^/'^). We compute 
characteristic properties of the photoemission spectrum 
at large finite A, such as the evolution of the height and 
position of the photoemission peak as a function of A. 
We thereby quantify the interpolation between the pho- 
toemission spectrum from strongly-coupled plasma and 
that from weakly-coupled plasma. 

Photoemission rate and spectral function. M — 4: SYM 
theory is a supersymmetric gauge theory with gluons, 
fermions and scalars all in the adjoint representation of 
SU{N), and a (global) i?-symmetry group S'[/(4). To 
model electromagnetism in this theory, one of the U{1) 
subgroups of the i?-symmetry group is gauged with cou- 
pling e. The Lagrangian can then be written as Q: 



L = LsYM + eJ^'^A''-\T\ 



(1) 



where Lsym is the Lagrangian oi N — 4 SYM theory, 
and the interactions internal to this Lagrangian are gov- 
erned by the 't Hooft coupling A — gyM^- The field 
is the photon (introduced by hand), with kinetic term , 
and is the electromagnetic current. The number of 
photons produced by a thermally equilibrated plasma per 
unit time per unit volume is given by 



dT. 
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where k is the 3-momentum of the (on-shell) photon, 
nb(/e) = Ijie^l^ — 1), T is the temperature, and aem = 
e^/47r. This holds to all orders in A, and to leading order 
in e. The quantity (k) is the light-like spectral density 
of the plasma defined via the retarded electromagnetic 
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current correlator Rn,u as 



-2Imi?^i,, where 



J3) 

Q{t) is the step function and K = (k, k), with k = \k\. 

We wish to compute R^^, using holography, but we 
do not have the holographic dual of the theory defined 
by Eq.dl]); we do, however, have the holographic dual of 
pure A/" = 4 SYM. The crucial point is that, to leading 
order in e, the electromagnetic current J^™ = J„, where 
Ji, is purely the i?-symmetry current associated with the 
U{1) subgroup. Therefore, the two-point function of the 
electromagnetic current can be replaced by the two-point 
function of the i?-symmetry current, computed entirely 
within the Af = A SYM itself Q . Our aim in this paper 
is to compute the retarded correlator of the i?-symmetry 
currents oi J\f — A SYM at strong 't Hooft coupling, 
retrieve the spectral function Xp,u, and insert the latter 
into Eq.([2]) to obtain the photoemission rate. Given that 
we wish to keep A strong yet finite, we must therefore 
work in the supergravity dual to M — 4 SYM plasma 
extended by the addition of finite a' corrections. 

Type IIB string theory setup at 0{a'^). The type IIB 
supergravity action corrected to 0{a'^) is Sub = Sjjg + 
Sfjg, where 
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where </> is the dilaton, the five-form, and i?io is the 
curvature. The leading 't Hooft coupling corrections are 
contained in Sfjg , and given schematically by 0, Q 
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where 7 = iC(3)(a7i?')^ = iC(3)A-3/2, with i?^ = 
gyj^Na'^ and C, is the Riemann Zeta function. C is the 
lOD Weyl tensor, and Tabcdef is defined as 



bcdef 
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abcnin def 



0171+ 771+ mn 

abfnin dec 



where the RHS is antisymmetrized in [a, 6, c] and [d, e, /] 
and symmetrized with respect to interchange of abc O 
def 0. Also, F+ = i(l + *)F5. The operators in Eq.® 
are dimension-eight operators obtained by various inde- 
pendent contractions of C and T that can be found in 5] . 
The background solution to this action is the corrected 
metric GmatQ, with f{u) = 1 — u^, 

" (5)'^ {-f{u)K\u)dt'+dx') 

P^{u)du^ + R^L''{u)dnl, (6) 



ds-" 



4w2/(w) 

where dfti is the line element on the S^, and 



K{u) 
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where the exponents are give by the expressions 
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(7) 



The extremality parameter is ro — uTR^ /{I + 2657/16), 
where T is the physical equilibrium temperature of the 
plasma. The boundary of the AdS space is at m = 0, and 
the horizon of the black hole at u = 1 . In addition, both 
-Fs and the dilaton have non-trivial background solutions, 
but their explicit forms are of no consequence in this 
work, as we explain shortly. The crucial point is that the 
tensor T is zero for the background solution Q . 

The vector perturbation. We now pursue the usual 
recipe involved in all holographic computations: firstly, 
perturb the supergravity background along the directions 
-0 which are dual to the field theory operators J whose 
correlation functions we are interested in, and plug the 
perturbed background into Sub- This yields the action 
S{'tp) of the perturbation ■)/'• Then, solve the equations of 
motion (EOM) of S{tj}) subject to -0 = "00 on the bound- 
ary of the space u = 0, and evaluate the on-shell action 
for these solutions, giving Z{ipo), the generating func- 
tional of correlation functions of the operators J^. Differ- 
entiating Z{iIjo) twice with respect to V'o yields < J. J >, 
and we are done. The details of this prescription in real 
time are described in Q. For the present case the per- 
turbation field ■0 dual to the i?-symmetry currents of 
the 4D theory is the vector perturbation of the grav- 
itational background obtained as a solution of the EOM 
of Sii B ■ The vector perturbation perturbs the metric 
and the Fc, solution, yielding 



ds'^ = gmn dx"" dx"- 



d^ij + iiKd(t>i + -^Af.dx'') 
V3 



(8) 



where gmn = Gmn for to, n S [0,4] and fii are the direc- 
tion cosines for the sphere, and 



R^L{u)- 

~7T~ 



J2 dp^r ^d(j}A^ *F2 , (9) 



where F2 — dA is the Abelian field strength of and 
e is a deformation of the volume form of the metric of 
the AdSBH- The Hodge duals * and * are taken with 
respect to the lOD metric and bH-AdSBH metric, re- 
spectively. Notice that the 7 — > hmits of Eq.dH [9]) 
are known exactly, see for instance [Sj. Both of Eq.(|8l 
[9]) are correct to linear order in 7. We will insert these 
Ansdtze into Sub = S'jjg + Sfjg below, obtaining an 
effective Lagrangian for which is at most quadratic 
in An. Two important points must be stated to this 



3 



end, dictated by the fact that we work strictly to hn- 
ear order in 7. For insertion into Sjj^, we require the 
Ansdtze to hnear order in 7. However, Sjjg only con- 
tains quadratic powers of F5 and, therefore, the e part of 
the Ansatz cannot contribute to the quadratic effec- 
tive Lagrangian of . On the other hand, Sfj^ contains 
operators which are higher than quadratic in F^, , so here 
e can contribute to the quadratic Lagrangian for A^, but 
the crucial point is that Sf^^ contains an explicit factor 
of 7 already, and so we only require e to zero order in 
7. This is of course nothing but the volume form on the 
AdS space Q . The incredibly simplifying upshot of these 
observations is that we do not require the 0(7) terms in e 
for our computation. This is why we do not care to state 
the explicit form of e. All we need to know is lim^^os- 
The same observations can be made for the contribution 
to the effective Lagrangian of A^ of operators containing 
the dilaton (and in fact any other field which is trivial in 
the zero order background supergravity solution) . These 
statements make the following work possible. 

The effective Lagrangian of the vector perturbation. 
Without loss of generality, we may set the photon mo- 
menta to (fc, 0, 0, k). In order to study the photoemission 
rate we only need to consider the transverse fluctuation 
Ax{t,x,u). Inserting the Ansdtze of Eq. dH [9|) into Sub, 
and integrating out the S^, we obtain a complicated La- 
grangian for Ax{t, X, u). This action can be coverted into 
a more useful form by use of the following field redefi- 
nitions: write = Ak{u)/{^y f{u)[l +^p{u)]), where 
p{u) is a polynomial in u beginning at O(m^) and Ak{u) 
is the Fourier transform of Ax{t, x, u). This takes us into 
the Schrodinger basis, such that the action is given by 
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(10) 



where /:*(u) = is the EOM *"(u) = y(u)^'(u), where 
the Schrodinger-like potential is given by 
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-^20984827/2 - 4752055U'* + 1838319u^ + q^u 
X (-16470 + 245442^2 + 1011173u'')] ) , (11) 



and q = fc/(27rT). The boundary term can be simplified 
to $ = ^''(m)^'(m) . We solve the Schrodinger equation 
analytically for the region k <^T, the so-called hydrody- 
namic regime of the plasma, and the high-energy regime 
k ^ T. In the intermediate momentum region we re- 
sort to a numerical solution of the EOM. Once we have 
the solution of the Schrodinger equation, the trace of the 
spectral function is given by: 



N'T' 
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The EOM of the vector perturbation. We solve the 
Schrodinger equation using perturbation theory. Write 
'^{u) = 5*0 (u) -I- 7\l/i(u), and insert into the Schrodinger 
equation, separating the powers of 7. The equation for 
^'o(m) is given by *^,'(m) = {-~f-''{u){l + q'u))-^o[u), 
and solves to give 
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*o(u) = (1 - u)"5(i+»9)(i + u)-i(i+ 

( {l + i)q {l + i)q ^ . l-u 
2Fi\l , , 1 - iq, — — 



(13) 



The equation for ^i(u) is solved numerically (if neces- 
sary). The trace of the spectral function x^(fc) is 



N'T' 
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which is exact to linear order in 7. We note that xjl(^) 
at A — )■ 00 is known [3], so our task here is to compute 
the 't Hooft coupling corrections to that result. 

Asymptotics of the spectral function, x'llik) can be 
evaluated analytically for low- and high-momentum, and 
numerically for the remaining momentum domain. We 
do not discuss the details of the computations, referring 
the reader to Q, and we simply display the results: 



l]\[2rp2 



( (1 + i^7) 9 + 0{q^) q « 1 
^ 1 4^St (1 + 57)9^/^ + 0(1) q»l 



with ( 1 
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7j coming from the factors of tq in Eq. ([T 



(14) 

The coefficient of q in the low-g regime of Eq. p^ means 
that the electrical conductivity of the strongly-coupled 
plasma is enhanced by a factor (l -I- 1^|^7) due to the fi- 
nite A corrections @ . This is as expected from the pertur- 
bative computations in Q: the weakly-coupled plasma 
has a larger mean-free-path per collision, allowing more 
efhcient diffusion of electric charge, and hence a higher 
electrical conductivity. The high-g region however, poses 
a question that we (as of yet) cannot answer: the authors 
of claim that the spectral function at weak coupling 
should go like q^^' in the ultraviolet (UV). Given that 
the spectral function at A — > 00 goes like q'^"^, in that 
regime one would have expected our result in Ea. (IT4l) to 
display some smooth interpolation between q^^' and q'^^. 
We do not obtain such an interpolation, finding instead 
that the finite coupling corrections do not change the 
g-dependence in the UV. Moreover, we find an enhance- 
ment by a factor (1 -I- 57) in that regime (see also fio}). 
The fact that the leading q behavior is unchanged by the 
corrections could have been seen from the Schrodinger- 
like potential in Eq. (fTT|) : the only g-dependence is q', 
identically to the A — >■ cx) case. Terms like g^, which 
could have changed the high— g functional dependence of 
Xftil), vanish. We shall revisit this point below. 

The photoemission spectrum. We plug the obtained 
X^(/c) into Eq.([2]) to give the photoemission spectrum. 
We display the results in figure 1. Clearly, the corrected 
result at strong coupling approaches the weakly-coupled 
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FIG. 1. The photoemission rate dV^Jdk in units of 
O.OlaemiV^r^ for different values of A, as a function of k/T. 
Solid, dashed, and small-dashed lines correspond to decreas- 
ing values of A oo, 75 and 50, respectively. The dotted 
line to the extreme left is the weak-coupling result at A = 0.5 
taken from Q. 



result (taken from reference [3|]). Firstly, the corrected 
curves exhibit a steeper tangent at the origin, due to the 
enhancement of the electric conductivity by the factor 
(1 + i^7) in Eq.dUl). Secondly, the peak of the pho- 
toemission is enhanced by the corrections, and the mo- 
mentum of maximal emission shifts towards the infrared 
(IR) , taking the corrected curves closer to the weakly cou- 
pled result. Simple numerical analysis on the light-like 
spectral function yields that the maximal rate is given by 



~dk 



0.0156695 1 



1115.3- 
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7j+0(7'), 
. . 

in units of a^mN^T^ , where we have made explicit the 
factor —265/87 coming from the overall normalization of 
the action. The position of the peak k„iax is 



1.48469 (1 - 842.4257) T + 0(7^ 



(16) 



This quantity can analytically be shown to be indepen- 
dent of the overall normalization of the action, making 
it an excellent candidate for comparing disparate gauge 
theories. One more quantity which is of interest is the 
total number of photons emitted, given by the area under 
the curves in figure 1. This is enhanced by a factor 



^total(7) 
^total(O) 



1 



461.941- 
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7 + 0(7'), (17) 



due to the fact that the weakly-coupled theory dominates 
in the IR, where Bose-suppression (due to n{,(fc)) is small. 

We finally make two comments about the behavior of 
the photoemission rate for high q. Firstly, there is a (A- 
independent) crossover point around k/T ^ 2.92, where 
the corrected curves dip below the A — 00 result. This 



is expected from the weak-coupling computations of [3|. 
What is surprising, as we mentioned above, is that the 
asymptotic values of the A-corrected curves for large k/T 
are given by (1 -I- 57) times the infinite coupling result, 
as in Eq. ([T^ (note that the domain of figure 1 does not 
extend to cover this asymptotic behaviour). This means 
that the finite-A corrections enhance the photoemission 
rate in the deep UV regime, contrary to the expectations 
of . Obviously, we are not guaranteed that the weakly- 
coupled result should be approached by strongly-coupled 
corrections computed in perturbation theory, especially 
not for a situation where the functional dependence on 
momenta is expected to be different, so we are not unduly 
concerned by this apparent discrepancy. It would be very 
revealing to understand these cross-over points, as well as 
their scaling with A. An important extension of our work 
would be to determine if the universality found in [llj for 
the energy-momentum spectral functions operates for the 
i?-symmetry current spectral function computed here. 
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